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I.  INTRODUCTION 


This  report  describes  the  results  of  an  in-house  program  conducted  within 
the  Air  Force  Rocket  Propulsion  Laboratory  (AFRPL).  The  objective  of  the 
program  was  to  develop  a  quick-response  rocket  motor  preliminary  structural 
analysis  capability  within  the  AFRPL.  Task  I  of  the  project  involved 
implementation  of  preliminary  structural  analysis  equations  on  the  HP  98I5A 
calculator.  Task  II  provided  the  means  to  calculate  an  effective  modulus  input 
for  thermal  transient  loading  based  on  linear  viscoelasticity,  and  provided  for 
graphical  and  numerical  output  utilizing  a  revision  of  the  linear  ramp 
approximation  program  used  in  the  original  task.  Task  III  provided  a  graphical 
and  manual  data  input  capability  to  allow  more  flexibility  and  Improved  accuracy 
In  the  overall  analysis. 


2.  SUMMARY 

2.1  "Quick-Look"  Motor  Structural  Analysis  Software.  The  preliminary 
structural  analysis  methods  given  In  The  JANNAF  Structural  Integrity  Handbook 
(Ref.  I)  and  the  Galclt  101  Technical  Report  (Ref.  2)  were  Implemented  In  an 
interactive,  user-oriented  program  for  the  HP  98I5A  calculator.  These 
equations  were  used  to  calculate  the  maximum  values  of  the  stresses  and  strains 
In  the  motor  bore  and  case-grain  bondline.  This  type  of  analysis  is  usually 
sufficient  to  determine  whether  or  not  a  given  grain  configuration  is  feoslblle 
and  If  a  more  detailed  analysis  is  required.  Comparing  the  preliminary  analysis 
to  more  complex  detailed  analyses  gives  the  analyst  a  reasonable  Idea  of  where 
the  results  should  be  falling. 

2.2  Thermal  Transient  Viscoelastic  Response  Software.  Analyses  Involving  a 
thermal  load  on  a  rocket  motor  require  an  effective  modulus  for  the  propellant 
at  any  given  time  (or  temperature).  The  approach  used  approximates  a  time- 
temperature  history  by  a  modified  power  law  that  results  In  a  tractable  form  of 
the  hereditary  Integral  to  approximate  the  stress  response.  Linear  strain  and 
temperature  coupling  were  assumed.  The  method  used  to  approximate  the  stress 
response  Involved  representing  the  strain  (as  a  function  of  reduced  time)  as  a 
series  of  linear  ramp  histories.  The  stress  response  was  calculated  for  each  ramp 
and  combined  with  each  proceeding  response  to  obtain  the  net  response  required. 
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2.3  Graphical-Manual  Input-Output  Facilitating  Software.  The  fact  that  much 
of  the  Input  data  required  for  the  above  analyses  is  available  only  In  graph  form 
established  the  need  for  these  programs.  To  complicate  matters,  many  of  these 
graphs  are  distorted  In  one  way  or  another  (photographs  taken  at  an  angle, 
photocopies  from  books  or  folded  originals,  to  name  a  few);  therefore,  a 
capability  in  the  graphical  Input  program  was  added  to  correct  for  these 
distortions.  The  other  Input  program  allows  for  manual  Input  of  data  (l.e.,  direct 
from  keyboard  of  calculator).  The  two  output  programs  allow  for  storage  of  data 
in  selected  registers  and  for  graphical  output  In  presentable  form. 

t 

2.4  Finite  Element  Analysis  Comparison.  Results  from  a  simple,  slow 
cooldown  problem  were  obtained  from  a  finite  element  analysis  using  the 
TEXGAP  2D  computer  code  and  the  "Quick-Look"  analysis  for  comparison  to 
verify  the  usefulness  and  accuracy  of  the  "Quick-Look"  analysis  and  supporting 
software. 

3.  DISCUSSION 

3. 1  "Quick  Look"  Motor  Structural  Analysis  Software.  The  JANNAF  Structural 
Integrity  Handbook1  contains  preliminary  structural  analysis  equations  that  can 
be  very  useful  for  a  quick  stress-strain  analysis  of  a  solid  propellant  rocket 
motor.  However,  use  of  these  equations  can  Involve  lengthy  hand  calculations. 

To  provide  a  usable  "Quick  Look"  analysis  tool,  It  was  decided  to  Implement  the 
equations  Into  a  set  of  software  for  use  on  the  HP  98t  5A  desk  calculator.  The 
equations  were  checked  for  consistency  and  some  necessary  alterations  were 
made  before  Implementa  i.  The  equations  were  also  modified  to  Include 
variable  propellant  and  ca»  no'  irial  Poisson's  ratio*.  1 

The  equations  varied  depending  on  the  type  of  motor  configuration  and  * 

analysis  performedi  therefore,  separate  program  files  were  created  for  each  type 
of  analysis.  Program  files  were  created  for  both  finite  motor  correction  factor 
determination  and  bore  stress  concentration  factor  determination  for  different 
bore  configurations.  These  program  files  are  used  as  routines  and  subroutine*  In 
a  main  program.  The  exact  equations  used  for  the  different  problem*  can  be 
found  In  Appendix  A.  The  aerodynamic  heating  analysis  and  the  thermal  load 
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analysis  use  the  same  equations,  except  the  temperatures  used  to  calculate  the 
shrinkages  are  different.  Specifically,  for  the  thermal  load  case, 
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and  for  the  aerodynamic  heating  case, 


8 


8 


c 


a 


(T  -  T,) 

^Tcase  '  V 


where  a  and  ac  are  coefficient*  of  linear  expansion  for  the  propellant  and  case, 
respectively}  TSF  Is  the  stress-free  temperature  for  the  motor}  T|  is  the  end 
condition  temperature  for  the  thermal  load  problem  and  the  beginning  propellant 
temperature  for  the  aerodynamic  heating  problem}  T  is  the  final  propellant 
temperature  in  the  aerodynamic  heating  problem  (T  is  usually  taken  to  be  equal 
to  T|)j  and  Tcase  Is  the  final  motor  cose  temperature. 


Figure  I  shows  the  temperature  histories  for  the  thermal  load  and 
aerodynamic  heating  problems  for  temperature  Input  comparison.  It  can  be  seen 
that  the  value  of  T|  is  the  same  for  both  analyses  for  a  given  problem.  Although 
the  propellant  temperature  normally  changes  somewhat  In  the  real  situation,  the 
value  T  is  taken  as  equal  to  the  value  T |  unless  the  user  specifies  a  new  value 
for  T, 


Figure  I,  Temperature  History  for  Thermal  Load  and  Aerodynamic  Heating  Problems 
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The  basic  equations  are  based  on  the  assumption  that  the  rocket  motor  is 
of  infinite  length.  This  may  be  a  valid  assumption  in  cases  where  the  motor  is 
long  and  small  in  diameter,  or  if  the  ends  are  bonded  to  the  case,  because  the 
end  effects  may  be  negligible.  However,  for  finite-length  motors  (i.e.,  stubby 
and/or  free  ends),  correction  factors  are  needed.  These  correction  factors 
("P  ,  Py,  and  Pj)  may  be  determined  from  graphs  found  in  The  JANNAF 
Structural  Integrity  Handbook'.  However,  accurate  visual  interpolation  from 
these  graphs  is  difficult,  because  the  graphs  are  three  dimensional  (two 
independent  variables).  Using  a  numerical  interpolation  method  has  two  distinct 
advantages.  Primarily,  fairly  good  accuracy  is  achieved--the  only  error  that 
occurs  Is  due  to  the  fact  that  the  curves  cannot  be  perfectly  described  by  the 
interpolating  functions.  The  approximations  are  acceptable  for  engineering 
accuracy  and  all  large  errors  are  corrected  (Figs.  2  -  4).  Second,  due  to  internal 
storage  and  numerical  Interpolation,  the  user  need  not  look  up  values,  which 
saves  a  considerable  amount  of  time.  Specifically,  the  method  used  is  two 
dimensional  Lagrange  polynomial  interpolation  that  approximates  the  surfaces 
created  by  the  graphs  with  a  polynomial  of  two  Independent  variables.  This 
methc  squires  data  from  the  graphs  to  be  set  up  In  an  Internal  matrix.  Enough 
data  had  to  be  used  to  adequately  represent  the  graph  (usually  a  7  x  7  matrix  was 
sufficient).  This  data  required  a  significant  amount  of  memory  space;  therefore, 
each  graph  used  a  separate  data  file  on  tape.  The  software  manipulates  the  data 
files  in  order  to  conserve  memory  space. 
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Figure  2.  Correction  Factor  -  Pr 


The  basic  equations  are  set  up  for  motors  with  circular  ports,  and  do  not 
include  stress  and  strain  concentrations  occurring  in  the  tips  of  slotted  ports.  To 
incorporate  these  effects,  equations  were  implemented  to  calculate  the  stress- 
strain  concentration  factors  for  varying  slot  tip  configurations.  The  equations 
used  wore  originally  derived  from  elastic  models  undergoing  infinitesimal  strains. 
However,  in  regions  of  high  s*ress  concentration  in  viscoelastic  material,  large 
deformations  can  be  expected.  These  large  deformations  tend  to  reduce  the 
stress  concentration,  because  the  material  deforms  Into  a  more  ideal  structural 
shape.  Therefore,  the  concentration  factors  obtained  tend  to  be  on  the 
conservative  side. 

The  acceleration  and  vibration  analyses  were  not  Implemented  at  tne  time 
of  this  writing.  The  main  program  lists  these  anlayses  as  choices  and  the  file 
exists,  but  the  programs  are  still  under  development.  The  thermal  load, 
aerodynamic  heating,  and  Ignition  pressure  analyses  were  Implemented  first 
because  of  similarities,  and  due  to  the  fact  that  these  analyses  were  required 
most  often.  The  acceleration  and  vibration  analyses  Involve  dynamic  loads  that 
require  a  different  procedure. 

After  the  software  was  tested  far  consistency  with  sample  problems,  It  was 
decided  to  allow  for  user  Input  of  the  Poisson's  ratio  for  both  case  and 
propellant.  The  equations  found  in  The  JANNAF  Structural  Integrity  Handbook  I 
assumed  a  Poisson's  ratio  of  0.50  for  the  propellant  and  0.33  for  the  case. 
Derivations  were  not  included;  therefore,  another  source  was  raquired.  The 
Galcit  101  Technical  Report^  contains  equations  with  variable  Poisson's  ratios. 
These  equations  were  used  to  obtain  equations  of  the  form  found  In  The  JANNAF 
Structural  Integrity  Handbook!. 

The  capabilities  of  the  "Quick  Look"  Structural  Analysis  software  are 
outlined  belowi 

a.  Stress-Strain  Analysis  Selection.  Any  or  all  cf  the  following  analyses 
may  be  selected: 

Thermal  Soak 
Aerodynamic  Heating 
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Acceleration  Effects  -  under  development 
Ignition  Pressure 

Vibration  Effects  -  under  development 

b.  Various  Simple  Grain  Configurations.  The  user  may  select  from  the 
following  configurations! 

Circular  Port 
Slotted  Port 

The  following  slot  tip  configurations  may  be  selected!  semicircular, 
ellptical,  or  flleted  flat-bottomed. 

c.  Finite-Length  Correction.  For  values  of  a/b,  correction  factors  are 
found  by  Lagrange  polynomial  Interpolation  of  Internal  graphs.  If  no  finite- 
length  correction  Is  necessary  the  correction  factors  are  assigned  to  the  value  of 
one  (I). 

d .  Superposition.  After  all  the  analyses  have  been  performed,  the  final 
program  file  sums  all  corresponding  stresses  and  strains  to  obtain  a  net  result. 

3.2  Thermal  Transient  Viscoelastic  Response  Software.  A  stress-strain 
analysis  of  a  solid  propellant  rocket  motor  undergoing  a  thermal  load  requires 
that  an  effective  modulus  of  elasticity  be  known  for  any  given  time  (or 
temperature).  To  obtain  an  effective  modulus  for  a  propellant,  some  knowledge 
Is  required  of  the  viscoelastic  processes  that  occur  In  the  propellant  under  a 
thermal  load.  One  of  the  complexities  Involved  In  the  analysis  of  a  viscoelastic 
material  Is  the  change  of  viscosity  with  change  of  temperature.  With  a  change 
In  temperature  of  a  material  with  linear  viscoelastic  properties,  a  time  shift 
occurs  In  the  response.  To  obtain  a  continuous  net  response,  the  individual 
curves  (plotted  on  log-log  coordinates)  are  shifted  so  they  fall  on  a  single  master 
curve.  The  time-temperature  shift  factors  are  determined  at  specific 
temperatures  by  applying  a  step  load  to  a  specimen  at  a  constant  temperature. 
The  response  is  then  recorded  and  plotted  on  log-log  coordinates.  This  process  Is 
repeated  at  the  specified  temperatures  to  obtain  each  response.The  time- 
temperature  shift  factors  can  then  be  determined  visually  or  by  numerical 
methods.  A  modified  power  law  was  used  to  represent  the  time-temperature 


shift  factors  as  a  function  of  temperature.  The  "reduced  time"  is  then 
determined  using  the  "Moreland-Lee  Hypothesis,"  the  representation  of  the  time- 
temperature  shift  factors,  and  the  temperature  history.  The  "Moreland-Lee 
Hypothesis"  states  that  for  a  nonconstant  temperature  history,  the  reduced  time 
can  be  represented  by  an  Integral, 
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Appendix  Q  contains  an  analytical  equation  for  the  reduced  time  that  was 
derived  for  use  in  the  software  along  with  the  derivation. 

Using  the  assumption  the  strain  and  temperature  are  linearly  coupled,  the 
strain  history  In  terms  of  reduced  time  can  be  computed.  A  method  was  then 
required  to  determine  the  stress  response  utilizing  the  hereditary  Integral. 
Polynomial  representation  and  piecewise  linear  ramp  approximation  were  the 
approaches  considered.  A  polynomial  representation  of  the  strain  history  would 
have  allowed  a  direct  solution  of  the  stress  response.  However,  the 
representation  was  too  complex  and  Inaccurate  and,  therefore,  was  not  used. 
Representation  of  the  strain  history  by  piecewise  linear  ramps  proved  viable  and 
accurate  and  was,  consequently,  the  method  used.  Basically,  the  method 
involves  breaking  up  the  strain  history  Into  a  series  of  linear  ramp  histories.  Tho 
stress  response  for  each  ramp  Is  computed  and  added  to  each  proceeding 
response  In  order  to  obtain  the  approximate  net  stress  response.  The  output 
values  of  stress  are  calculated  at  offset  values  of  reduced  tlmo 
0.e.,  _i+kCJ(£-|‘£i_i )  where  0<_k  <_1],  which  results  in  a  more 

accurate  representation  of  the  stress  response.  kaIs  found  by  trial  and  error  (a 
typical  value  for  Is  .833),  and  should  be  fairly  non-sensitive  if  a  large  number 
of  ramps  are  used. 

To  evaluate  the  hereditary  Integral,  the  relaxation  modulus  must  be  known 
for  any  given  reduced  time.  The  relaxation  modulus  data  Is  normally  given  as 
the  net  response  obtained  from  the  propellant,  relaxation  tests.  A  modified 
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power  law  fit  was  used  once  again  to  represent  the  data  in  a  form  the  calculator 

could  easily  handle.  If  the  temperature  normalized  form  of  the  relaxation 

modulus  data  Is  given  [i.e.,  E  ■  (^  )  Er3,  then  the  computed  stresses  must  be 

corrected  by  multiplying  the  stress  values  by  corresponding  factors  of  (X  ) . 

'o 

From  the  values  of  stress  and  corresponding  values  of  strain,  an  effective 
modulus  (a/e)  may  be  computed  at  each  output  point.  Values  of  real  time  can 
also  be  computed  by  using  the  same  equations  that  were  used  to  compute 
reduced  time. 

The  capabilities  of  the  thermal  transient  viscoelastic  response  software  are 
outlined  belows 

a  *  Representation  of  Relaxation  Modulus  Data.  Modified  power  law  fit 
by  a  least  squares  Iterative  method. 

b.  Representation  of  Time-Temperature  Shift  Factor  Data.  Modified 
power  law  fit  through  three  points. 

c.  Representation  of  Temperature  Rate.  Modified  power  law  fit 
through  three  points.  Exact  representation  of  constant  temperature  rate  (linear 
temperature  history). 

d.  Plot  or  Storage  for  Ramps  Approximation.  Plot  or  storage  of  strain 
history  (In  terms  of  reduced  time).  Plot  nontime  also  allows  visual 
representation  of  ramp  approximation  for  use  as  a  check. 

e .  Allows  for  Correction  of  Stress. 

f .  Choice  of  Plot  or  Table  of  Data.  Allows  choice  of  parameters  to  be 
plotted  or  printed  on  hard  copy.  Data  is  not  stored.  Parameters  are  as  foilowsi 

Real  Time,  t 
Reduced  Time,  5 
Strain,  * 

Stress,  u 

Effective  Modulus,  E(£) 
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3.3  Graphical-Manual  Input-Output  Facilitating  Software.  Data  needed  for 
analysis  may  come  in  several  different  forms  that  require  special  attention. 
Data  in  tabulated  form  usually  can  be  handled  easily  by  direct  input.  However, 
data  in  graphical  form  can  be  difficult  to  use,  particularly  if  the  graph  contains 
distortions  due  to  photocopy  errors,  etc.  The  HP  9862A  plotter  has  a  unique 
feature  that  allows  the  user  to  digitize  graphs.  Basically,  the  user  places  the 
plotter  in  the  digitizer  mode,  moves  the  pen  (pointer)  over  the  point  of  Interest, 
and  then  exits  the  digitizer  mode.  The  coordlnaies  of  the  point  are  entered  into 
the  calculator  automatically.  However,  the  digitizer  assumes  the  graph  is 
undistorted.  Therefore,  for  distorted  graphs,  corrections  need  to  be  made.  The 
method  used  to  correct  distorted  graphs  was  Lagrange  shape  functions  (Ref.  3). 
This  technique  maps  an  undistorted  grid  defined  by  the  user  into  the  distorted 
grid,  thereby  creating  an  internal  distorted  grid,  which  when  digitized,  will 
produce  correct  data.  This  data  may  then  be  replotted  in  a  presentable  form 
using  the  plotting  routine  (Fig.  5),  or  stored  in  user  selected  registers  using  the 
data  manipulation  routine.  The  plotting  routine  utilizes  a  cubic  spline  subroutine 
(Ref.  4)  to  accurately  fit  sparse  data.  The  plotting  routine  will  accept  data  from 
either  the  manual  Input  program  or  the  digitizing  program. 


0,0  0,4  0,8  1,2  1.6 


TIME,  MIN  X  10-3 
Figure  S.  Plot  F.xampie. 
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The  capabilities  of  these  programs  follow. 

a .  Data  Input 

(1)  Manual  Input.  Direct  Input  of  data  Into  calculator 
through  keyboard.  Used  when  tabulated  numerical  data  Is  available. 

(2)  Digitize/Correct.  Input  of  data  by  way  of  the  HP  9862A 
Plotter/Digitizer.  Used  when  data  Is  In  graph  form  only.  Distorted  graphs  (l.e., 
perspective-linear  and  warped-curved  distortions)  are  corrected  by  using  a 
technique  that  utilizes  Lagrange  shape  functions. 

b.  Data  Output 

(1 )  Allows  user  to  define  the  memory  area  In  which  to  store 
the  inputted  data  for  use  elsewhere.  Data  may  be  stored  on  tape. 

(2)  Plot.  Allows  user  to  plot  data.  Allows  for  any  size  graph 
(to  plotter  llm!ts)i  axis  headings,  and  graph  titles. 

3.4  Finite  Element  Analysis  Comparison 

To  verify  the  validity  and  usefulness  of  the  "Quick  Look"  structural 
analysts  software,  a  simple,  slow  cooldown  problem  was  run  on  this  software  and 
on  the  TEXGAP  2D  finite  element  computer  code  (Ref.  5).  The  problem 
description  and  results  can  be  found  In  Appendix  C. 

Several  plots  were  made  to  show  the  accuracy  of  the  representations  of  the 
time-temperature  shift  factors  and  the  relaxation  modulus  rate  (Figs.  6  and  7). 
Figure  6  also  shows  the  control  dote  points.  These  points  were  the  data  used  In 
the  program  to  represent  the  time-temperature  shit  i  factors.  As  can  be  seen, 
these  points  He  on  the  curve,  and  the  mid  and  final  points  are  true  data  points. 
The  first  point  Is  not  an  actual  data  po!nt|  It  was  chosen  so  tne  resulting  curve 
would  besr  represent  the  actual  data. 


Figure  6.  Time-Temperature  Shift  Factor  vs  Temperature. 


Figure  7.  Propellant  Relaxation  Modulus  vs  Reduced  Time. 

Plots  were  also  made  of  the  final  results  from  the  "Quick  Look"  structural 
analysis  (Figs.  8  -  10).  These  results  In  numerical  form  were  compared  to  the 
results  found  with  the  TEXGAP  2D  code.  The  "Quick  Look"  (structural  analysis 
produces  stresses  and  strains  theoretically  located  at  the  bore  surface  and 
bondllne  In  the  axial  midplane  of  the  motor,  while  the  TEXGAP  2D  code 
produces  stresses  and  strains  throughout  the  motor  case  and  grain  that  Is  defined 
by  an  element  grid.  Therefore,  the  values  obtained  from  the  TEXGAP  analysis 
were  taken  from  elements  whose  sides  were  located  on  the  axial  midplane  of  the 
motor.  The  average  percent  difference  (using  the  "Quick  Look"  structural 
analysis  results  as  the  base  values)  for  the  axlsymmetrlc  problem  was  around 
4  percent  to  6  percent  (Table  I).  The  difference  In  the  results  may  be  attributed 
to  many  reasons.  Namely,  these  analyses  are  determined  by  two  totally 
different  approaches.  The  "Quick  Look"  structural  analysis  equations  are  a 
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direct  method  for  obtaining  preliminary  results.  The  TEXGAP  2D  code  is  an 
analysis  designed  to  give  the  analyst  a  detailed  view  of  the  stress  state  in  a 
motor  undergoing  a  given  load,  but  it  is  not  necessarily  going  to  give  highly 
accurate  results  at  specific  points  within  the  motor,  since  the  accuracy  depends 
on  the  fineness  of  the  finite  element  grid.  The  "Quick  Look"  structural  analysis 
gives  theoretically  exact  results,  but  errors  can  result  in  the  determination  of 
the  correction  factors.  As  mentioned  previously,  the  correction  factors  are 
determined  through  an  Interpolation  technique  and  are,  therefore,  approximate. 
Because  of  this  uncertainty  In  the  determination  of  the  correction  factors,  a 
second  analysis  was  done  using  the  same  motor  parameters  as  before,  but  with  an 
infinite  length.  The  average  percent  difference  was  around  1.3  percent  to 
5  percent  (Table  2).  The  largest  percent  difference  occurred  In  the  values  of  the 
hoop  strain  (5  percent).  These  percent  differences  are  within  acceptable 
engineering  standards. 


Figure  8.  Hoop  Strain  vs  Real  Time 
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Table  1.  Comparison  of  “Quick  Look"  and  TEX6AP  Structural  Analysis  Results  (Finite-Length  Motor) 
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APPENDIX  A 


Equations  Used  In  the  "Quick  Look"  Structural 


Analysis  Software 


Thermal  Load 


If  ^  0,5(  ob  -  ob. 

If  p  >  OiSj  Ujj  *  Ojj 

•  2 

but  with  an  equivalent  Cp 


<f »«  ■  t^2  -  <S»J 


where  subacript  A  denote*  calculation  of 

A  2 

(jp)  using  a  radius  "a"  of  a  circular  port 
of  equivalent  cross  sectional  area  to  the 
slotted  port. 


°T  *  ^b 

where  H  is  a  stross  concentration  factor 
for  slotted  ports. 
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Thermal  Load 


Thermal  Load 


°T  "  Hob 

O'  -  (o')v  i-  B6 
z  t 

ca  "  Vr  C1  *  +  ^ 
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Pressure  Load 


PROBLEM  DESCRIPTION 


D  FINITE -LENGTH, 
CIRCULAR-PORT  GRAIN 
WITH  INTERNAL  PRESSURE 
LOADING -ENDS  BONDED 


Sene  bquetions  es.  for  Problem  I 


Pressure  Load 


-  (i  -  Tr>r 

-  (l  -  fyp 


r  v  "^‘wHwet-v  *  umtu 


Pressure  Load 


PROBLEM  DESCRIPTION 

U 'LONG  SLOTTED  PORT  GRAIN 
WITH  INTERNAL  PRESSURE 
LOADING  OU  FINITE-LENGTH 
GRAIN  WITH  ENDS  BONDED 


f  i  °*5.  °b  “  ab* 

If  £  >  0.5,  ab.  p.  cb 

o  2 

but  with  an  equivalent  Og-) 


(|)2  -  0.2  [$2  - 


where  subscript  A  denotes  calculation  of 

a  2 

(gO  using  a  radius  "a"  of  a  circular  port 
of  equivalent  cross  sectional  area  to  the 
slotted  port. 


oT  -  H(P  +  ob)  -  P 

where  H  is  a  stress  concentration  factor 
for  the  slotted  ports. 


Pressure  Load 


PROBLEM  DESCRIPTION 


10UG  SLOTTED -PORT  GRAIN 
WITH  INTERNAL  PRESSURE 
LOADING  OR  FINITE-LENGTH 
GRAIN  WITH  ENDS  BONDED 


Pressure  Load 


PROBLEM  DESCRIPTION 


FIN  ITE -LENGTH, SLOTTED -PORT 
GRAIN  WITH  INTERNAL 
PRESSURE  LOADING-ENDS 

FREE 


•  IF  ENDS  ARE  DOME-SHAPED 
ADD  1/2  OF  DOME  HEIGHTS 
TO  L. 


o'  -  H(P  +  o')  -  P 

o'  -  HV(P  +  O')  -  P 

C'  -  F  [|(p  ♦  ob)Cl  -  V2)] 


APPENDIX  B 


Derivations  of  Equations  Used  In  the 
Thermal  Transient  Viscoelastic  Software 
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APPENDIX  C 

JANNAF  linear  viscoelastic  Round  Robin  ana  I /sis  problem 
description  and  comparison  results  from  the  "Quick  Look"  analysis 
and  the  TEXGAP  finite  element  analysis. 

The  description  contains  two  problems.  Problem  I,  only  the  slow 
cooldown  problem  was  analyzed}  therefore,  references  to  Problem  2 
may  be  Ignored. 
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JANNAF  STRUCTURES  AND  MECHANICAL 
BEHAVIOR  SUBCOMMITTEE  (S&MBS) 

JANNAF  EXPERIMENTAL  AND  ANALYTICAL  STRUCTURAL  METHODS  PANEL 
LINEAR  VISCOELASTIC  ROUND  ROBIN 

I.  BACKGROUND 

A.  Introduction 

.  Problem!  Several  methods  of  analyzing  solid  propellant  motor  stresses 
and  strains  based  on  linear  viscoelastic  theory  are  In  use  by  propulsion  contractors 
and  government  agencies.  These  Include  "quasl-elastlc"  methods  (In  which  elastic 
analyses  by  finite  element  codes  are  used  together  with  certain  mathematical 
procedures  to  provide  viscoelastic  response)  and  "linear  viscoelastic"  computer 
codes.  The  degree  to  which  these  methods  yield  the  same  results  when  applied  to 
the  same  problem  Is  unknown. 

B.  Oblectlve 

The  objective  of  this  "round  robin"  analysis  Is  to  collect  results  from 
different  analysts,  using  various  techniques,  on  the  same  quasl-statlc  linear 
viscoelastic  problem.  These  results  will  be  tabulated  and  reported  to  the 
JANNAF  S&MBS.  The  overall  goal  Is  to  determine  how  consistent  the  various 
analytical  procedures  are  -  i.e.,  do  they  all  yield  similar  results? 

C.  Ground  Rules 


Confidentiality  of  the  source  will  be  maintained  In  all  reporting  of 
results.  The  source  of  Individual  analysis  results  will  be  known  only  to  the 
reviewers  and  will  be  used  only  for  follow-up  questions  and  for  recognition  of 
participation.  Participants  may  elect  to  analyze  either  one  (or  both)  of  the  two 
problems  posed  below.  Each  participant  will  be  provided  a  copy  of  the  report.  By 
stating  Intent  to  participate,  each  participant  agrees  to  carry  out  the  analysis 
according  to  the  problem  statement,  to  provide  a  report  of  the  results,  and  to 
respond  to  subsequent  reasonable  queries  for  additional  Information.  The 
voluntary  nature  of  this  participation  Is  clearly  recognized,  however. 

40 


II.  PROBLEM  DESCRIPTION 


A.  This  round  robin  analysis  consists  of  two  motor  analysis  problems: 

Problem  No.  I.  A  "slow”  cooldown  (slow  enough  to  maintain  thermal 
equilibrium  throughout  the  propellant  grain). 

Problem  No.  2.  A  "rapid"  cooldown  problem,  In  which  thermal 
equilibrium  cannot  be  assumed. 

B.  Analysis  approaches  may  Include  your  normal  techniques  or  other 
methods.  The  propellant  shall  be  treated  as  a  "Thermorheologically  simple 
material."  Other  simplifications  to  propellant  material  behavior  may  be  applied  as 
well.  All  assumptions  and  simplifications  (e.g.,  numerical  integration  techniques) 
shall  be  reported. 

C.  The  problem  details  are  given  In  Tables  I  through  4.  Tobies  I  through  3 

provide  data  that  are  common  to  both  problems.  Table  4  gives  the 

temperature-time  "histories"  for  Problems  No.  I  and  2,  respectively. 

D.  The  results  shall  be  reported  on  the  "Report  of  Results"  forms  (Tables  4 
and  5).  A  narrative  report  is  also  encouraged. 

TABLE  I.  Geometry,  Constraints  and  Material  Properties. 
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Propellant  Properties 

E  (see  Propellant  Relaxation  Data)  (Table  2) 
v  *  0.499 
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TABLE  2.  Propellant  Relaxation  Data. 
Propellant  Stress  Relaxation  Modulus 

Reduced  time,  min  Ere.’  psl 


I0“3 

1,300 

I0-2 

I0_l 

900 

630 

10° 

460 

I01 

355 

I02 

290 

I03 

250 

I04 

220 

I05 

200 

I06 

ISO 

I07 

168 

I08 

150 
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Time-Temperature  Shift  Factor 


Temperature,  °F  Log  Ay 

-45  7.52 

-42  6.60 

-29  5.19 

-2  3.84 

19  2.20 

76  0.00 

118  -1.23 

152  -2.30 


TABLE  3.  Additional  Data. 


Steel  thermal  properties* 


p  «  0.282  10/In3 
c  *  0.127  BTU/ft-F 
k  *  266  BTU/hr  ft2-F 


Propellant  thermal  properties* 


p  *  0.0636  lb/ln3 
cp  =  0.298  BTU/ft-F 
k  a  0.275  BTU/hr  ft2-F 
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TABLE  4.  Report  of  Results  for  Slow  Cooldown 


Problem  No.  I.  Slow  Cooldown  (2®F/day). 

Elapsed  time,  Temperature,  Hoop  strain  ee, 
days  °F  '"/in 
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